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Abstract 


Bats are notorious reservoirs of genetically-diverse and high-profile pathogens, and are 
playing crucial roles in the emergence and re-emergence of viruses, both in human and 
in animals. In this report, we identified and characterized previously unknown and 
diverse genetic clusters of bat coronaviruses in the Atlantic Forest Biome, Brazil. These 
results highlight the virus richness of bats and their possible roles in the public health. 


Brazil harbours 15% of the world’s bat diversity and carries 178 distinct bats species 
(Nogueira MR, 2014). Out of these, a total of 113 species exist only in the Atlantic 
Forest Biome (AFB), which is the second largest rain forest of the South America, and 
is one of the unique regions with highest biodiversity in the world (Paglia AP, 2012). 


Bats are historically unique and widespread mammals playing essential roles in the 
emergence and re-emergence of viruses of both veterinary and public health 
importance. Although viruses of diverse genetic backgrounds can co-exist 
asymptomatically in bats, majority of these viruses are single stranded RNA viruses 
(Calisher et al., 2006). Coronaviruses (CoVs) are enveloped, positive-sense, single- 
stranded RNA viruses in the family Coronaviridae, and are usually associated with 
respiratory, enteric, hepatic and neurological pathologies of varying severity (Woo et 
al., 2012). Coronaviruses are classified into four genera; A/phacoronavirus (O-CoV) and 
Betacoronavirus (B-CoV) have been exclusively identified in mammals, whereas 
Gammacoronavirus and Deltacoronavirus are mainly detected in avian species (Woo et 
al., 2012; ICTV. 2015). Based on the genetic relatedness, the B-CoVs can further be 
subdivided into four clades: A to D (Drexler et al., 2014). 


All CoVs that can potentially infect human were originated from animal reservoirs, and 
four of such CoVs are believe to be transmitted to human through bats including a 
CoVs (229E and NL63) and highly pathogenic B-CoV (Severe Acute Respiratory 
Syndrome and Middle East Respiratory Syndrome) (Bolles et al., 2011; Chan et al., 
2015; Corman et al., 2013; Huynh et al., 2012). Recently, a number of novel bats CoVs 
have been identified, primarily from African, Asian and European bats (Calisher et al., 
2006; Chu et al., 2006; Drexler et al., 2014), as well as from South American countries 
including Costa Rica, Panama, Ecuador, Mexico and Brazil (Corman et al., 2013; Goes 
et al., 2013). Collectively, these studies indicate the co-existence of bats and viruses at 
the interface of viral evolution and bats ecology. 


A limited number of studies have been conducted in Brazil to map the nature and breath 
of bats in harbouring viral populations. In previous studies, a total of five distinct CoV 
lineages have been detected in just 10% of local bats (15 species), and most of these are 


belonging to a-CoV (Brandao et al., 2008; Corman et al., 2013; Goes et al., 2013). 
These attributes, and the existence of a large number of human beings (120 million) in 
the Atlantic Forest Biome, Brazil, clearly highlight the potential of bats in not only 
carrying zoonotic viruses but also in possible transmission of viruses to human beings. 


To ascertain the diversity of bats CoVs circulating in the Brazilian bats, a total of 401 
intestine tissues from 17 bat species were collected from 2010. to 2014 encompassing 
bats with distinct diet habits (Table 1, Figure 1A, B). Bats from urban (N=192) and 
rural (N=10) areas were received from 14 municipalities of the Northwestern State of 
Sao Paulo by the Rabies Laboratory of Universidade Estadual Paulista (UNESP), 
Aracatuba-SP, and Zoology Department at UNESP, Rio Claro-SP. Additionally, 99 and 
100 samples were collected from Iguacu National Park, and from two distinct disturbed 
landscape sites, respectively. All activities were authorized and approved by the Ethic 
Committee of Institute of Biomedical Research from University of Sao Paulo (56-18- 
03/2014) and bats species were identified based on the morphological characteristics 
including head-boby, forearm size, and dental arch as described by Vizotto & Tadei, 
1973; Gregorin & Taddei 2002; Miranda et al., 2011. 


Nucleic acids extracted from 30 mg of intestine tissue of individual bat using 
NucliSENS® easyMAG® were subjected to first strand cDNA synthesis with random 
hexamers and RT-PCR High-Capacity cDNA Archive Kit (Applied BioSystems). The 
cDNA preparations were screened by pancoronavirus Nested PCR assay targeting the 
RNA-dependent RNA polymerase (RdRp) gene (Chu et al., 2011). Positive samples 
(n=15) were sequenced and final sequences were submitted to the GenBank 
(KT717380-KT717394) (Supplementary table). A dataset, consisting of sequence 
generated in this study and all publically available bats CoVs, was phylogenetically 
analysed by neighbour-joining method in MEGA6 software using Kimura’s two- 
parameter correction and 10000 bootstrap values (Tamura et al., 2013). 


Coronaviruses were detected in 15 bat intestines samples from eight bat species with 
distinct diet habit, demonstrating a marked potential of CoVs distribution among bat 
species in AFB that harbours 9% of world’s bat diversity. All CoVs-positive bat species 
were geographical distributed in the Neotropics and in anthropogenic area apparently 
affected by fragmented forests, and bats-abundance urbanized areas. The detection of 
CoVs varied by bat species; undisturbed forest remnant (4%), fragmented forests (4.5%) 
and urban areas (3.1%) (Table 1). 


Phylodynamics analyses indicated the circulation of 13 a-CoV and 2 B-CoV in AFB 
bats (Figure 1C). This is of primary interest as majority of bat CoVs surveillance studies 
conducted in the New World bats detected only aCoVs (Corman et al., 2013; Drexler 
et al., 2014; Osborne et al., 2011). A higher resolution analysis indicated a distinct 
distribution and diversity of 0-CoV and B-CoV lineages (Figure 1D). 


a-CoV sequences obtained from bats of same genus presented high nucleotide sequence 
similarity (e.g. Artibeus, Glossophaga, Carollia, Molossus, Myotis and Sturnira) (Figure 
1D and Supplementary table), even with sequences detected in other studies from bats 
of geographically distant regions. This relation can be exemplified by the high similarity 
of the CoV RNA partial sequences (91,9%) detected in Carolla perspicillata species 
from Fenix-PR, Brazil (KT717385), and Fyzabad, Trinidad and Tobago (EU769557) 


(Carrington et al., 2008), located at a distance of at least 3,900 Km. Similar results were 
previously reported for a variety of bat CoVs and are taken as evidence of co-evolution 
of CoV genotypes and specific host genera (Drexler et al., 2014; Corman et al., 2013; 
Anthony et al., 2013). 


Moreover, two previously uncharacterized O-CoV from Myotis nigricans and Myotis 
riparius and one B-CoV cluster were identified in studied bats (Figure. 1C and D). CoV 
sequences from Myotis bats genera presented high nt sequence similarity (98.9%) and 
were most closely related to bat CoVs (KC110771) collected from Brazilian Molossus 
molossus and Tadarida brasiliensis. A distinct cluster of O-CoVs were also detected in 
Sturnira bat that grouped with a-CoVs lineage 1, a group with an evolutionary history of 
recombination and cross-species transmission between domestic and livestock animals, 
such as feline, canine and swine a-CoVs (Lorusso et al., 2008). This clade presents nt 
similarity between 73.9% and 74.4% with transmissible gastroenteritis virus 
(DQ811789) and feline coronavirus (AF124987) depending on the Sturnira lineage 
(Supplementary table). This branching pattern may possibly explain the common- 
ancestral origin of GCoVs lineage 1 species from bats and other animal species. 
However, extensive evolutionary studies on complete genome sequences of these 
isolates are required to provide information on the virus origin and divergence. 


Out of two B-CoVs, the virus from Eumops glaucinus clustered within MERS-CoV 
containing lineage C. The second B-CoV detected in Artibeus /ituratus shown highest 
similarity (96.4%) with CoV reported from Costa Rica (Corman et al., 2013) and lower 
similarity (69.3%) with MERS-CoV (Supplementary table). These sequences showed 
sequence similarities with B-CoV detected in Pipistre/lus bats that poses high homology 
with human MERS-CoV (JX869059). Notably to observer that the Fumops bat positive 
for CoV was found on an urban area and was predated by a domesticated cat. It was not 
possible to rule out any event of virus transmission to human, however, this highlight 
the potential of domesticated animals in virus transmission and disease dynamics at the 
virus-animal-human interface. This is of special importance due to the established role 
of animals in the transmission of viruses to human being (Johnson et al., 2015). Taken 
together, these results represent the first detection of lineage C B-CoV in South 
American bats. Despite of close relationship between lineage C B-CoVs in Asia, Africa 
and Europe, cumulative data indicate that this lineage showed a more diversified host 
family distribution in Americas in Mormoopidae, Phyllostomidae and Molossidae bat 
species (Anthony et al., 2013; Corman et al., 2013; Goes et al., 2013). 


Conclusions 


We present a great diversity of CoV genotypes and clusters in Brazilian bats, 
highlighting a biogeographic distribution of bats CoVs in the region. It is indispensable 
in future to investigate the evolutionary events in genetically diverse bats CoVs using 
complete genome sequences, and their possible transmission potentials to human being. 
Although it is not possible to calculate the risk of “spill over” events of Brazilian bats 
CoVs to humans, our results reinforce the need for expanded and continuing 
surveillance of CoVs in bat fauna, including those in the AFH regions of Brazil. 
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Figure legend 


Figure 1. A: Brazil’s map and sites for bats capture/collection. Bats were collected from 
17 sites (a-q) in the Atlantic Forest Biome from two adjacent Parana and Sao Paulo 
states, Brazil. B: Based on morphological characteristics, bats were classified into 
different families (positive families were bold and underlined). Phylogenetic tree was 
adapted from Simmons, N.B. 2005. C: Phylogenetic analysis of bat CoVs detected in 
the AFB, and reported previously around the world. Brazilian sequences reported in this 
study are marked with arrows. D: A higher resolution tree of the alpha and 
betacoronaviruses. Sequences studied in this report are marked with black filled square. 
Original map available in http://mapas.sosma.org.br from Atlas of the Atlantic Forest 
Remains (Fundagaéo SOS Mata Atlantica/Instituto Nacional de Pesquisas Espaciais- 
INPE). 
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Highlights 


1. We describe the circulation of CoV lineages, 13 a-CoV and two B-CoV in AFB bats 


2. Coronaviruses were detected in 15 bat samples from eight bat species (50% 
positivity) 


3. Acluster of a-CoVs were detected in Sturnira bat that grouped with a-CoVs lineage- 
1 


4. We identified two distinct B-CoV that clustered within MERS-CoV containing lineage 
C 


5. Results highlight the virus richness of bats and their possible roles in the public health 


